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INTRODUCTION 


PURPOSE 

This report summarizes the results of effort aimed toward the 
development of a rudimentary model that illustrates the relative mass 
requirements of food production components in a controlled ecology 
life support system (CELSS) based on regenerative concepts. The report 
is intended as a working paper which can provide a basis for further 
model development and analysis. 

The model and analytic results can be useful for developing an 
understanding of the mass requirements for food production in a CELSS 
and how these requirements compare with food resupply requirements. 

Such an understanding aids in making knowledgeable decisions about 
research investment options in regenerative life support. More 
importantly, the documentation of the model and results reveals gaps 
in knowledge and thereby provides guidance for improving the model and 
the analysis procedures. 

BACKGROUND 

The need for life support options based on regenerative methods 
becomes increasingly important from a logistics standpoint for long 
duration missions. This need has been recognized and discussed in 
several workshops and reports (1,2, 3, A, 5). 

Other studies, reports, and papers have proposed using scenario 
analyses to help evaluate research progress and to help program managers 
assess technology options for regenerative systems (6,7,8,9,10,11,12). 

The scenario analysis approach, although utilized in these prior studies, 
has heretofore not been documented in the detail necessary for compara- 
tive studies by other researchers. This report provides documentation 
of the initial model and illustrates its use and its sensitivity to 
changes in assumptions and parameter values. In addition to document- 
ing the overall scenario model, this report also updates the diet 
assumed for the 1977 summer study (e.g., see 9). 
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MODEL SCOPE AND STRUCTURE 

SCOPE 

Assumptions 

The model is a simple deterministic formulation that equates the 
mass required for food supply under two different scenarios: A) complete 

resupply of diet materials (baseline scenario) and B) part of the diet 
supplied by food grown in the space habitat by regenerative methods 
(alternative scenario). The model thus uses time as the dependent 
variable, permitting the calculation of a "breakeven time," which may 
be considered as the mission duration for which the total mass required 
for supplying dietary requirements under the baseline scenario equals 
the mass required under the alternative scenario. 

Several simplifying general assumptions are implicit in the 
resulting formulation: 

1. Mass requirements for the baseline :food supply scenario 
are linear with time. 

2. Total mass requirement for each scenario may be represented 
by a simple summation of component masses. Each component 
represents a particular necessary function and may further 
be subdivided into expressions representing subcomponents 
or macrocategories of design factors. 

3. Several functional requirements which may be substantially 
different in the two scenarios are included from the model: 

- energy collection, power generation, energy storage, 
and waste energy rejection; 

- human and food supply waste processing requirements; and 

- food and waste storage requirements. 

4. Food production by conventional biological growth (plant 
photosynthesis) is a viable candidate food production 
process. (See also reference 13.) 

Focus and Limitations 

The above assumptions indicate the model's focus on the food 
production function itself and not on components and functions which may 
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be heavily dependent on the choice of food supply. The model, at this 
stage of development, therefore is useful primarily for comparing con- 
ventional food supply alternatives which are judged not to have a sub- 
stantial difference in requirements for (or impacts on) energy and waste 
processing functions. The results of model calculations also can indicate 
orders of magnitudes of mass requirements on an absolute scale, providing 
a starting point for assessing alternative food supply scenarios and 
for determining the relative mass of food supply components. 


MODEL STRUCTURE 


Baseline Scenario (Resupply) 


The baseline scenario is characterized by full resupply of necessary 
diet materials. The diet requirements are often given on a person-day or 
person-year basis. By letting by the resupply mass required for food 
(g/person-year) , the food resupply mass for a population of N p for a 
period of T years is given by 

(M R )(N p )(T). (1) 

To this mass should be added the mass required for air revitalization. 
Engineering considerations (14) suggest that this mass can be modeled by 
an expression representing a fixed mass and plus a mass which is linearly 
dependent on the number of persons: 

M o + «ov><V- (2> 

The mass required for the resupply is thus given by the sum of (1) plus (2) 


Alternative Scenario 

The alternative scenario is characterized by recycling of nutrients, 
utilizing plants for growing food. The mass required for growing and 
processing the food may be represented by the sum of five mass components: 
the food-producing biomass (M^) , the mass of the harvesting equipment 
the food production system mass (system components for plant growth, 
excluding the plant biomass itself, Mpg) » the mass of the food processing 
equipment (M^) , and the mass of the water of transpiration (M^) 
required in the plant growth environment (including the atmospheric/vapor 
phase and the nutrient reservoir but excluding the water contained in the 
plant biomass). Each mass would be given in grams (or kilograms) unless 
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( 3 ) 


other dimensions are specified. The total mass for the alternative 
scenario is given by 

^ ^ + + ^FP + “wT* 

Calculations for each of the components are given in the following 
paragraphs. 

Food Producing Biomass . The total biomass required for the food 
production function is given by the sum of the biomass required for each 
of the individual foods in the diet. Letting be the biomass required 
to produce the daily dietary requirement of one person for food j , the 
total productive biomass is given by 




(4) 


is calculated from the dietary requirements which are met by the 
food produced by the plants. If the daily diet amount of food j is R . , 
this is equal to the production rate of fresh edible food (g/person-day) 
and thus the total biomass harvest rate, B . (g/person-day) is given by 

^oj = (5) 

where e . is the edible fraction (fresh). Assuming a simplified model 
of linear biomass growth over time, with t^ as the total growing period 
(time to harvest), then 


My - <v<V 2 >- (6) 

Harvesting Equipment . Food harvesting equipment for a CELSS has not 
been designed, and much of the harvesting functions may be fulfilled 
manually. Grain products (e.g., wheat) may be an exception, and this 
model assumes a single grain harvester (plot harvester) with a mass of 
800 kg (15). That is, the initial assumption for harvesting equipment 
requirements is: 

Mjj = 800. (7) 

Food Production System . The food production system consists of the 
components necessary to sustain productive plant growth in a controlled 
environment. These components include lighting, atmospheric control, 
plant support, and control components. Estimates for these components, 
based initially on a design for earth-based phytotron facilities, have 
been placed into three categories: (1) fixed mass components, whose 

mass is independent of the total growing area required (e.g., control 



components) , (2) an unsealed mass which is linearly related to the size 
of the growing area (e.g., lighting components), and (3) a scaled mass 
which increases with the size of the growing area by a power relation- 
ship (e.g., atmospheric control components). The last category uses the 
empirical engineering rule of 0.6 scaling, in which the mass requirement 
increases with the 0.6 power of the area. Such a scaling is consistent 
with components whose capacity would be related to volume rather than 
area, for example. (Appendix A gives the individual components of a 
ground-based controlled environment (phytotron-type) system, their mass, 
and the estimated masses for the three model categories.) 

The model calculation for the food production system mass (Mp s ) is 
therefore 


^PS ^PSC + ^ M a^ ^ + ^ M as^ ’ (8) 

where is the total fixed (control component) mass, A is the growing 

area (or total illuminated area) , M is the factor for the unsealed 
(linearly related to area) mass, and M is the scaled area factor. 

2 

The biomass growing area required per person (A , in m /person) is 

D 

calculated from the dietary requirements for the diet components and the 
plant productivities. For food j, the required area is given by 

— bj - . (9) 

where is the dietary requirement (required fresh edible production 

rate), e^ is the edible fraction (fresh) of the total biomass harvested 

daily, and £. is the biomass productivity, or total biomass growth rate 

2 ^ 

(g/m -day) . 

Food Processing Equipment . Food processing equipment mass is 
estimated from data on current food preparation system technology and 
assumed engineering developments (16) . The mass assumed for this 
equipment is 


^ = K(N p ) X , (10) 

where the factor K and the power x are determined by fitting a curve to 
the data for populations of 10 and 100 persons, and these parameters are 
302.7 and 0.415, respectively, for Mpp in kilograms. 

Water Reservoir Mass . The final component of mass in the alternative 
scenario is the water of transpiration, (kg). This mass includes the 
atmospheric (vapor phase) moisture in the plant growth chamber and the 
mass of the nutrient reservoir, but excludes the water in the plant 
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biomass. Plants differ in the amount of water they transpire daily, but 
the amount can be modeled by a linear relationship to the dry biomass of 
the plant. Consequently, in order to compare different plants and 
different diets, the model incorporates a parameter which represents 
the number of days of transpired water assumed to be required in the 
system for the growing biomass. For plant j, the water of transpiration 
mass, is given by 


Vi ' doHVjXV'V- (11> 

where is the dry biomass (kg/person) of plant j, 10 is the empirical 

factor representing the ratio of water mass transpired (per day) to the 
plant's dry biomass, is the population, and is the design factor 
representing the number of days of transpiration water required by the 
system. 


L 
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EXAMPLE CALCULATIONS 


DIET AND PLANT GROWTH ASSUMPTIONS 
Diet Basis 

The diet requirements assumed for both the baseline and alternative 
scenarios are based on a "thrifty food plan." This diet is characterized 
by less consumption of animal tissue and somewhat more grain consumption 
than actual consumption patterns of food stamp recipients. The diet 
assumed for the model is the thrifty plan with small variations due to 
the assumed percentage of waste and to minor diet variations. Table 1 
compares the actual food consumption pattern, the thrifty diet, and the 
diet assumed for the model. Note that the last column in the table also 
is equivalent to the required production rate of edible fresh food. 

Plant Growth Rates 

Plant growth data are taken from the literature on crop yield and 
plant productivity studies. The other model parameters related to ratios 
of fresh and dry weights, edible and total biomass, etc., are taken from 
the same literature. Table 2 presents these and the related calculated 
values for the major food plants in the assumed diet. 


STORAGE AND RESUPPLY ASSUMPTIONS 

For the baseline scenario, food must be stored or resupplied. The 
mass of this stored or resupplied food is assumed to be linearly related 
to the edible dry biomass, with the small residual moisture content and 
packaging mass being 53% of the dry mass. The daily resupply mass for 
food j is therefore 




= (1.53)EB 


dj’ 


( 12 ) 


where EB^ is the edible biomass, dry, of food j. 

The mass required for air revitalization, a function assumed to be 

performed adequately by the plants in the alternative scenario, is given 

by expression 2 for the baseline scenario. From reference 14, the total 

mass for air revitalization has a small (90.4 kg) fixed component and a 

larger variable component. The parameters for expression 2 are 

M = 90.4 and M = 415.8. 
o ov 
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Table 1. Food Quantities - Comparison of 
Consumption Pattern, Thrifty Diet, and Model Assumptions' 1 ' 

(fresh weight, grams/person-day) 


Food 


Consumption 
Pattern ^ 


Thrifty 
Food Plan 


Diet Assumed for 
Model Calculations 

<V 


A. Milk, cheese, ice cream 

548.4 


404.5 


399.8 

B. Meat , poultry , fish 

209.2 


196.5 


196.5 

C. Eggs 

41.1 


34.3 


30.2 

D. Dry beans, peas, nuts 

16.2 


28.5 

Beans 

16.2 




Peanut Butter 

11.7 

E. Potatoes 

84.3 


131 


131 

F. Dark green, deep 

24.3 


25.3 

Leaf Cabbage 

13.6 

yellow vegetables 




Carrots 

11.7 

G. Citrus fruit , tomatoes 

110.3 


116.7 

Tomatoes^ 

116.7 

H. Other vegetables , fruit 

230.2 


239.3 

Green Beans 

44.8 





Head Cabbage 

12.3 





Lettuce 

20.1 





Me Ions 5 

253 





Peas 

9.7 

I. Grain products 


Cereal 

57.7 

Wheat Equiv. 

352 



Flour 

59.7 





Bread 

148.5 




Other 

Bakery 





Products 

86.3 



J. Fats, oils 

38.9 


61.6 


61.6 

K. Sugar, sweets 

55.1 


55.8 


55.8 

L. Accessories Not Given 

Not Given 

80.4 


Notes : 

^"Adapted from "The Thrifty Food Plan" (17) and "Report to J. Spurlock" 
(BSSG) from Marcus Karel (18) . 


2 

Based on National Survey of Food Stamp and Food Distribution Program 
Recipients, November, 1973 (Average 4-person Household); from "The 
Thrifty Food Plan" (17). 

3 

Based on 20-5A year-old male nutritional requirements; amounts for 
some groups allow for approximately 5% discard and waste (e.g., 
egg shells) . 


Tomatoes replace oranges on equal mass basis. 


5 


Melons replace apples and bananas on 2. 48 g per gram basis. 
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Table 2. Plant Productivity Parameters for Major Foods 


Column 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Reference 

* 

Parameter 

R 

e 

t H 

TB, 

d 

P 

EB d 

EB f 

e B 

B 

— o 

4 

w 

e> 


Food 







(6*4) 

(1-7) 

1* (4x7) 

l-(5*6) 


Dry Beans 

16.2 

47 

43.1 

204.3 

21.1 

24.3 

.119 

136.1 

.67 

.132 

19 

Peanut Butter 

12.3 

110 

49.5 

355.0 

8.2 

9.4 

.026 

473.1 

1.33 

.126 

20 

Leaf, Head Cabbage 

25.9 

30 

10.4 

180.8 

9.9 

172.1 

.952 

27.2 

.15 

.942 

21 

Carrots 

11.7 

80 

43.6 

315.6 

21.3 

154.2 

.489 

23.9 

.08 

.862 

21 

Tomatoes 

116.7 

215 

7.9 

113.5 

6.2 

95.0 

.84 

138.9 

1.22 

.93 

22 

Potatoes 

131.0 

120 

20.2 

134.6 

13.7 

97.9 

.727 

180.2 

1.34 

.86 

23 

Green Beans 

44.8 

60 

108.9 

920.6 

26.3 

305.2 

.332 

134.9 

.15 

.914 

24 

Lettuce 

20.1 

28 

11.6 

221.1 

8.5 

161.3 

.730 

27.5 

.12 

.947 

25 

Melons 

253.0 

107 

32.9 

396.5 

19.9 

298.9 

.754 

335.5 

.85 

.933 

26 

Peas 

9.7 

50 

15.6 

99.9 

0.6 

3.7 

.037 

262.2 

2.62 

.838 

27 

Wheat 

352.2 

196 

148.4 

505.3 

58.5 

67.2 

.133 

2648.1 

5.24 

.129 

24 

Totals 

993.6 


486.1 





4387.6 

13.77 




R = daily requirement (g-person-day , fresh); t = time to harvest (days); TB, = total biomass, dry; 

e H a _2 -1 
P = productivity; EB, = edible biomass, dry; EB- = edible biomass, fresh (all g-m-day ) ; e„ = edible 

d r -1 -1 2r -1 

fraction, fresh; B^ = biomass harvest rate (g-person-day ) ; Ag = growing area required (m-person ); 

w = fraction moisture, edible, 
e 



BREAKEVEN CALCULATION PROCEDURE AND RESULTS 

The calculation of times (years) for which the baseline and 
alternative scenario mass requirements are equal initially was per- 
formed using two programs for a programmable hand calculator (HP-67). 
These programs, described in Appendixes C and D, permit rapid evaluation 
of detailed scenarios and the evaluation of the model's sensitivity to 
changes in single parameter values. Appendix B summarizes the nomen- 
clature conventions for both programs. 

The first program (Appendix C) calculates the plant productivity, 
growing area requirements, and the other calculated plant parameters 
from the fresh and dry total and edible biomass productivity factors, 
the dietary requirements, and the growing period (time to harvest) for 
individual cultivars. 

The second program uses the results of the first and design 
parameters (population and assumed parameter values for the growing 
environment and food processing component masses) to calculate the 
mission duration for which the mass requirements of the two scenarios 
are equal. This program is a preliminary, incomplete version of the 
model outlined in the text above. Specifically, the program in 
Appendix D does not include expression (2), the mass assumed to be 
required for air revitalization components in the baseline scenario. 
Therefore, the program is particularly useful for calculating breakeven 
times for individual food items/cultivars (scenario comparisons in 
which air revitalization is required both for the baseline case and 
an alternative that involves limited on-board food production) . 

Table 3 presents the results of calculations for the major plant- 
produced foods. Note that wheat has the earliest breakeven time, 
followed by beans (dry and green), melons, and potatoes, with peas 
having the longest breakeven time. 

For the alternative scenario in which virtually all plant-derived 
food is grown, with the biomass providing both food production and air 
revitalization functions, the breakeven time is 12.4 years, as shown 
in Table 4. This table also shows the results of modifying the initial 
values of the parameters to generate additional alternative scenarios 
for comparison with the resupply case. 
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Table 3. Breakeven Analysis Results - Comparison of Resupply and Regenerative Scenarios 


Assumptions: Resupply mass = 1.53 x dry mass of food required. 

Initial alternative scenario parameter values: 

"p - 1° M fp ■ 0 Wpsc ' 15 - 88 M a s ' 122 




T = 1 
r 

“H = ° M a 

= 83.2 





(Baseline) 


(Alternative) 




Food 

Resupply Mass 

Biomass Holdup 

Growing Env. Equip. 

^WT 

T 

BE 



(kg-person-yr^) 

(kg-person^) 

Mass (M PS )(kg) 

(kg) 

(yrs) 

A. 

Dry Beans 

7.87 

3.22 

955.3 

69.9 

13.1 

B. 

Peanut Butter 

6.00 

25.6 

1677 

355.5 

34.3 

C. 

Cabbage 

.84 

.407 

361.8 

2.4 

43.4 

D. 

Carrots 

.90 

1.01 

189.2 

13.9 

22.7 

E. 

Tomatoes 

4.56 

14.9 

1578 

109 

37.3 

F. 

Potatoes 

10.24 

10.8 

1710 

162.3 

18.4 

G. 

Green Beans 

2.15 

4.14 

296.3 

46.6 

, 16.1 

H. 

Lettuce 

.60 

.402 

266.8 

2.1 

44.9 

I. 

Melons (for Apples) 

4.00 

7.64 

578.5 

63.4 

16.2 


(for Bananas) 

2.33 

4.46 

381.0 

37.0 

18.3 


(for Fresh Fruit) 

3.13 

5.94 

475.1 

49.3 

16.9 

J. 

Peas 

.88 

6.54 

3061 

102.1 

360 

K. 

Grain (Wheat) 

171.3 

259.5 

5688 

7629 

7.9 



Table 4. Breakeven Analysis, All Plants Combined 


Initial V alues 

Mpsc “ 15 * 9 ( k S) 


Hj. r = 214.8 (kg person-yr* - ) 

Np = 10 (persons) 

T = 1 
r 

M, = 787* 
fp 

Hg = 344.6 (kg person'*') 

Mg^ = 86.4 (kg person* - ) 

Mg = 800 (kg) 


M a = 83.2 (kg m 2 ) 

M as = 122 (kg m 2 ) 

2 -1 
Ag = 13.8 (m-person ) 

M = 90.4 
o 

M = 415.8 
ov 

= 0.415 in equation (10). 


Other values given in Table 2. 

* 

Corresponding to K = 302.7 and x 


Results 

Scenario Breakeven Time ( T ER » Years) 


A. 

Initial values, above 

12.4 

B. 

As 

A, 

with M„ = M =0 
o ov 

12.4 

C. 

As 

A, 

with peas resupplied 

9 

D. 

As 

c. 

with T r - 1.1 

8.5 

E. 

As 

c. 

with peanut butter resupplied 

9.7 

F. 

As 

c. 

with NP = 100 

8.5 

G. 

As 

E, 

with N = 100 
P 

7.7 


12 





DISCUSSION 


ISSUES FOR FURTHER STUDY 

The results of applying the model thus far indicate several issues 
which should be considered further. These issues are outlined in the 
following paragraphs. 

Plant Productivity Measures 

The area of the growth chamber has a relatively large impact on 
the overall CELSS mass requirements. Consequently, plant productivity, 
as measured by the yield of edible (and digestible) biomass per unit 
area, is an important consideration. 

Another contributor to total mass of the regenerative system is 
the standing biomass, or "biomass holdup." Plants which have low biomass 
holdup for a given edible production are thus desirable. (Note that 
time to harvest, or maturity period, is not directly, or solely, an 
adequate measure of productivity.) 

Functional Components 

The results indicate that the components associated with the growth 
chamber (scaled and unsealed components) are significant contributors 
to the total mass. Reductions in the mass required for these components 
(e.g., lighting) should be possible. 

CONCLUSIONS 

The model, although limited in scope, provides a useful starting 
point for analyzing alternative diet and plant growth scenarios. Further 
refinement is underway to improve the model's utility and to facilitate 
its application by other researchers. 
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Appendix A 

CONTROLLED ENVIRONMENT WEIGHT ESTIMATES 


System Component 


Initial 1 
Per .98m 2 


Estimated Weights (lbs) 


(Net) 

Fixed 


Optimistic^ 


(Net) 
Per ,98m 2 


(For .6 scaling) 
Per .98m^ 


1. Panels 612 

2. Fan Coil Assembly 150 

3. Refrigeration 

Cond. Unit 185 

4. Lighting 150 

5. Control Console 76 

6A. Humidification 33 

b. Dehumidification 205 

7A. Exhaust & Make up 

Air Blower 76 

B. CO 2 Analyzer 

8. Plant Support 

Structure 150 

9. Misc. Components 160 

(Wireways-50; 

Insulation-10; 

Piping-50; Vibration 
Links-5; VJires-30; 
Valves, etc. -15) 


25 


10 


(Include in 6B) 


(Include 

w/heat rejection) 
150* 


33 

200 


30 


30 


,0 

Total 

(lb) 

1797 

35 

180 

263 



(kg) 

812.12 

15.88 

81.65 

119.30 



(kg-m 2 ) 

827.6 

NA 

83.21 

121.57 


David Raper, personal communication (estimates :from Environmental Specialties, 
Raleigh, North Carolina). 

2 

Estimates by Mike Modell, Jack Spurlock, Dave Raper, and Bob Mason (79/06/26). 

•k 

Reduction should be possible. 
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Appendix B 


' NOMENCLATURE AND CONVENTIONS 

Conventions 

Subscript j refers to type of crop in cultivar. 
Total biomass defined as edible and inedible. 
Underbar: per person basis. 


Parameter or Variable 

-2 -1 

Growing area (m -person ) 

Total biomass harvest rate (g person-day^") 

Edible fraction, fresh 

Mass requirement for growing area, 
unsealed components (kg-m~2) 

Mass requirement for growing area, 
scaled components (kg-m~2) 

Dry biomass holdup (g-person^") 

Food processing equipment mass for 10 persons (kg) 

Mass of food harvesting equipment (kg) 

Total growing area equipment mass (kg) 

Mass of control equipment for food production module (kg) 

Additional mass required for alternative scenario 
above that required for baseline scenario 

Mass required for resupply/storage of food j 
(kg-person - l-yr - l) 

Mass of (transpiration) water required for' adequate 
humidity and nutrient reservoir 


Symbol 

'BJ 

M 

a 

M 

as 

M, 

fp 

"H 

^PS 

^sc 

Hj-p 

^TRj 


18 



Parameter or Variable (cent.) 


Symbol 


Population 

, -2 , - 1 . 

Total growth rate of fresh biomass (g m-day ; 

Production rate of fresh edible food (g personiday'*') 
(Set equal to daily requirement) 

Harvest waste, fresh (g person-day^") 

Breakeven time (yrs) ; mission duration for which mass 
requirements of baseline and alternative scenarios 
are equal 

Total time to harvest (days) 

Number of days of transpiration supply assumed 
Fraction of water in total biomass 
Fraction of water in edible portion 
Fraction of water in inedible portion 
Harvest waste water (g personiday^) 


- N 

P 


P. 


R . 


R . 

-wj 





w . 
ej 


- w . 
WJ 


W . 
wj 



Appendix C 

HP-67 PROGRAM FOR CULTIVAR CALCULATIONS 


Storage 


0 R . 

-ej 


1 w . 
ej 


2 w 


wj 


3 e Bj 


4 % 


6 JL • 

7 (TB dj ) 

8 < E V 

9 (EB ) 


Program Instructions 


1. Key in TB, . 

dj 


2 . 

3. 

4. 

5. 

6 . 

7. 


" " TB.. 

fj 

" " EB, , 

d I 

" " EB.. 

fj 

II II 

Hj* 

" " R . , 

-ej 

Read out : 


, press | ENT 


ENT 


1) 

R . 
~ej 

7) 

w . 
ej 

2) 

w . 
ej 

8) 


3) 

w . 
WJ 

9) 

B . 

-oj 

4) 

e Bj 

10) 


5) 

C Hj 

ID 

R . 

6) 

P. 

3 

12) 

W . 

-VJ 
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AGRICULTURAL DATA INPUT CONVERSIONS 


Input Data Arrays 


TB . : Total biomass, 

dj 

TB^ : Total biomass, 

EB : Edible biomass 

dj 

EB,..: Edible biomass 

fj 


dry (g m 2 day 

-2 -1 
fresh (g m day ) 

, dry (g nf 2 day 1 ) 

-2 -1 
, fresh (g m day ) 


Derived Variable Arrays 


w e j : Fraction of water in edible portion 

w Fraction of water in inedible portion 

wj 

e_ . : Edible fraction, fresh 

Bj 

-2 -1 

P.: Total growth rate of fresh biomass (g m day ) 


Sequence of Calculations 


w 


ej 


EB £ . - EB. . 

_Jj iL = i _ 


EB 


fj 


1 - EB. ./EB,. 
dj fj 


w . 
wj 


(TB fj - EB fj ) - (IB d1 - EB,.) - EB^ 

-(TB fj -EB f .) 


<TB £3 - 



EB 

TB 


11 

fj 


P. = TB. . 
J fj 
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ALGORITHM FOR AGRICULTURAL PRODUCTION 


Design Variables 


N 


Population 


R . : Production rate of fresh edible, food (g person ^ day 


Data Arrays 


w 


w 


w J 


'B'j 

C H3 


p. 

~3 


Fraction of water in edible, portion 
Fraction of water in inedible portion 
Edible fraction, fresh 
Total time to harvest (days) 

-2 -1 

Total growth rate of fresh biomass (g m day ) 


Derived Variables 


U Bj 

“Bj 

B . 
-°3 

% 

R . 
-wj 

W . 
“WJ 


Fraction of water in total biomass 
Growing area (m^ person 

Total biomass harvest rate (g person ^ day 
Total biomass holdup (g person 
Harvest waste, fresh (g person ^ day 
Harvest waste water (g person ^ day 


Sequence of Calculations 

1. w_ . = w . e. . + w . (l-e_.) 
Bj ej Bj wj B j 


2. A n . = R ./ P. e . 

-03 ~3 Bj 

3. B . = R ,/e D . 

-oj — ej Bj 

4. M_ . = B . t .12 

-Bj -oj Hj 

5. R . = B . (1-e . ) 
-^3 -oj Bj 6 7 

6. W . = R . w . 

— wj — wj wj 
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CULTIVAR CALCULATIONS PROGRAM LISTING 


1. 

* LBL A 

28. 

* 


55. 

- X 

- 

<v 

2. 

STO 9 

29. 

STO 3 


56. 

RCL 

0 


3. 

1 

30. 

R/S 


57. 

RCL 

5 


4. 

STO 8 

31. 

* LBL 

B 

58. 

T 



5. 


' 32. 

STO 0 


59. 

RCL 

3 


6. 

STO 5 

33. 

- X - 

(R .) 
— e i 

60. 

f 



7. 

+ 

34. 

+ 


61. 

- X 

- 

<V 

8. 

STO 7 

35. 

STO 4 


62. 

RCL 

0 


9. 

1 

36. 

RCL 1 


63. 

RCL 

3 


10. 

RCL 8 

37. 

- X - 

(w . ) 
ej 

64. 

i 



11. 

RCL 9 

38. 

RCL 2 


65. 

- X 

- 

( V 

12. 

V 

39. 

- X - 

( V 

66 . 

STO 

6 


13. 

- 

40. 

RCL 3 


67. 

RCL 

4 


14. 

STO 1 

41. 

- X - 


68. 

X 



15. 

RCL 7 

42. 

RCL 4 


69. 

2 



16. 

RCL 8 

43. 

- X - 

C V 

70. 

T 



17. 

- 

44. 

RCL 5 


71. 

- X 

- 

<v 

18. 

RCL 5 

45. 

- X - 

<v 

72. 

1 



19. 

RCL 9 

46. 

RCL 1 


73. 

RCL 

3 


20. 

- 

47. 

RCL 3 


74. 

- 



21. 

* 

48. 

X 


75. 

RCL 

6 


22. 

CHS 

49. 

1 


76. 

X 



23. 

1 

50. 

RCL 3 


77. 

- X 

- 

(R •) 
-wj 

24. 

+ 

51. 

- 


78. 

RCL 

2 


25. 

STO 2 

52. 

RCL 2 


79. 

X 



26. 

RCL 9 

53. 

X 


80. 

- X 

- 

(W .) 

— WJ 

27. 

RCL 5 

54. 

+ 


81. 

R/S 
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Appendix D 

HP-67 PROGRAM FOR BREAKEVEN CALCULATION 


Storage 


Register Parameter 
I N 


9 

8 
. 7 


M 


fp 


^R 


Register Parameter 

6 npsc 


5 

4 

3 


4 

4d 

4 


Register Parameter 

2 4 


i 

0 


M 


as 


M 


Operating Instructions 
1. Key in N„ 


Press ENT 


2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 


9. 

10. " " M. 

11 . 

— o 

12. Read out: 


4r 

T 

r 

% 
4 * 
4d ' 
4 
4sc 

M 


10 


-3 


10 


-3 


as 


m 

ENT | 


BE 


m . 

| ENT | 

II 

[I] 

Breakeven time (years) 


Notes 


(kg) 


,-3 


-1 


(Mg *10 ) = kg-person 

-3 -1 

(Mg^*10 ) = kg-person 


M tp Total PCELSS mass requirement (above resupply 
component mass) 


M 


wt 


Water mass - transpiration reservoir 

Mpg Total growing environment equipment mass 

To run again changing only a few parameter values: 

key in new value, press STO X , where X is from above table. 
Repeat for each changed value 


Press (GT^ Q (0] [|| |5| |[7s 
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CELSS BREAKEVEN ANALYSIS PROGRAM LISTING 


r 


i) 


* 


1. 

* LBL 

A 

29. 

ENT 

57. 

1 

2. 

STO 9 


30. 

RCL 0 

58. 

5 

3. 

4 


3i . 

X 

59. 

Y* 

A. 

STO 8 


32. 

X Y 
■*- 

60. 

RCL 9 

5. 

4 


33. 

• 

61. 

X 

6. 

STO 7 


3A . 

6 

62. 

STO C 

7. 

4 


35. 

Y X 

63. 

RCL 3 

8. 

STO I 


36. 

RCL 1 

63a . 
63b. 

RCL I 

w 

A 

9. 

R/S 


37. 

X 

6A . 

+ 

10. 

* LBL 

B 

38. 

+ 

65. 

RCL 5 

11. 

STO 6 


39. 

RCL 6 

66. 

+ 

12. 

4 


AO. 

+ 

67. 

RCL A 

13. 

STO 5 


Al. 

STO A 

68. 

+ 

1A. 

4 


A2. 

RCL I 

69. 

RCL B 

15. 

STO A 


A3. 

RCL 8 

70. 

+ 

16. 

4 


AA. 

X 

71. 

STO D 

17. 

STO 3 


A5. 

RCL A 

72. 

RCL 7 

18. 

R/S 


A6. 

X 

73. 

RCL I 

19. 

* LBL 

C 

A7. 

1 

7A . 

X 

20. 

STO 2 


A8. 

0 

75. 


21. 

4 


A9 . 

X 

76. 

- X - 

22. 

STO 1 


50. 

STO B 

77. 

RCL D 

23. 

4 


51. 

RCL I 

78. 

- X - 

2A . 

STO 0 


52. 

1 

79. 

RCL B 

25. 

RCL 1 


53. 

0 

80. 

- X - 

26. 

RCL 2 


5A . 


81. 

RCL A 

27. 

X 


55. 

• 

82. 

- X - 

28. 

ENT 


56. 

A 

83. 

RTN 


(Hr P ) 


(Mt P ) 
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Calculations 



i. 

"fss - M as <4V 6 



2. 

"psu ‘ *Wp 



3. 

= npss + ^su + 

^sc 

(STO A) 

A. 

M wt " 10 * *Bd * \ 

. T 

r 

(STO B) 

5. 

M = M &)* A15 

fp ”fp 4o' 


(STO C) 

6. 

h b = V4 



7. 

=m b + m h + m fp 

+ “ps + M wt 

(STO D) 


T BE ■ “tp * «TR • V 


8 . 



APPENDIX E 


FORTRAN Program for Breakeven Calculation 


Purpose: This program' is based on the PRELIMINARY SCENARIO 

ANALYSIS MODEL by Robert M. Mason. 

Author: Martha Sadler, New View, Sept 1980 


Environment: 

DEC VAX-11/780 VAX/VMS VAX-11 FORTRAN IV-PLUS 
Non-Standard Code: 

VAX-11 FORTRAN IV-PLUS extensions of ANS FORTRAN 1966: 

Data types: CHARACTER, LOGICAL* 1, 

Block IF logical structure: IF THEN, ELSE, ELSE IF, END IF 

statements 

END= and ERR= in READ or WRITE statements 

OPEN, CLOSE, DEFINE FILE file control specifications 


Commons Used: 


<name> 
REALl 
CHAR2 
INT3 
LOG 4 


<descr iption> 

Real variables 
Character variables 
Integer variables 
Logical variables 


Subroutines 
<name> 
BLOCK DATA 
PRODUC 
BREAK 
DESCRT 
MFREIN 
MDECDE 


Called: 

<descr iption> 

COMMON DATA 

This subroutine performs cultivar calculations. 
This subroutine performs breakeven calculations. 
Gives description of variables when '?' input. 
Modified version of FREEIN.FOR written by Walton 
Modified version of DECODE. FOR written by Walton 


Limitations 
Exit Points 

Constraints and Cautions 



VARIABLE DESCRIPTIONS 


TR NUMBER OF DAYS OF TRANSPIRATION 

MCONST RATIO OF PLANT THAT IS DRY WEIGHT 
TBDJ TOTAL DRY BIOMASS (G PER SQ M PER DAY) 

TBD TOTAL DRY BIOMASS (G PER SQ M PER DAY) 

TBFJ TOTAL FRESH BIOMASS (G PER SQ M PER DAY) 

TBF TOTAL FRESH BIOMASS (G PER SQ M PER DAY) 

EBDJ EDIBLE DRY BIOMASS (G PER SQ M PER DAY) 

EBD EDIBLE DRY BIOMASS {G PER SQ M PER DAY) 

EBFJ EDIBLE FRESH BIOMASS (G PER SQ M PER DAY) 

EBF EDIBLE FRESH BIOMASS (G PER SQ M PER DAY) 

WEJ FRACTION OF WATER IN EDIBLE PORTION 

MBD DRY BIOMASS HOLDUP (G PER PERSON) 

MPSC MASS OF CONTROL EQUIPMENT FOR FOOD PRODUCTION MODULE (KG) 
MA MASS REQUIREMENT FOR GROWING AREA, UNSCALED COMPONENTS 
(KG PER SQ M) 

MAS MASS REQUIREMENT FOR GROWING AREA, SCALED COMPONENTS 
(KG PER SQ M) 

MFP FOOD PROCESSING EQUIPMENT MASS FOR 10 PERSONS (KG) 

MH MASS OF FOOD HARVESTING EQUIPMENT (KG) 

WWJ FRACTION OF WATER IN INEDIBLE PORTION 
EBJ EDIBLE FRACTION FRESH 

PJ TOTAL GROWTH RATE OF FRESH BIOMASS (G PER SQ M PER DAY) 

WBJ FRACTION OF WATER IN TOTAL BIOMASS 

ABJ GROWING AREA (SQ M PER PERSON) 

REJ PRODUCTION RATE OF FRESH EDIBLE FOOD 
(G PER PERSON PER DAY) 

BOJ TOTAL BIOMASS HARVEST RATE (G PER PERSON PER DAY) 

THJ TOTAL TIME TO HARVEST (DAYS) 

TH TOTAL TIME TO HARVEST (DAYS) 

RWJ HARVEST WASTE, FRESH (G PER PERSON PER DAY) 

MBJ TOTAL BIOMASS HOLDUP (G PER PERSON) 

MB TOTAL BIOMASS HOLDUP (G PER PERSON) 

TOTMBJ BIOMASS HOLDUP FOR TOTAL POPULATION 

TOTMFP FOOD PROCESSING EQUIPMENT MASS FOR TOTAL POPULATION (KG) 
MTRJ MASS REQUIRED FOR RESUPPLY/STORAGE OF FOOD J 
(KG PER PERSON PER YR) 

MTP ADDITIONAL MASS REQUIRED FOR ALTERNATIVE SCENARIO ABOVE 
THAT REQUIRED FOR BASELINE SCENARIO 
MWT MASS OF WATER REQUIRED FOR ADEQUATE HUMIDITY AND 
NUTRIENT RESERVOIR 

MPS TOTAL GROWING AREA EQUIPMENT MASS (KG) 

NP POPULATION 

TBE BREAKEVEN TIME (YEARS) 
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FORTRAN PROGRAM 


C SEE 'PRELIMINARY SCENARIO ANALYSIS MODEL', CELSS WORKING PAPER, 

C WP754-1 BY ROBERT MASON, APRIL, 1980 

REAL DEFVAL (6 ,11) , NUMBER , AMOUNT , GENVAR ( 6 , 2 ) ,USRVAL {6 ,11) ,TR, 

X MCONST ,TBDJ-,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA,MAS ,MFP ,MH , 

XWWJ ,EBJ ,PJ ,WBJ , ABJ ,REJ ,BOJ ,THJ ,RWJ ,MBJ ,TOTMBJ ,TOTMFP , 

X MTRJ ,MTP ,MWT,MPS ,NP,TBE, ALLMBD, ALLMTR, ALLMBJ , ALLABJ 
CHARACTER* 13 FDNAME (11) , WHICH 
CHARACTER* 80 LINE 

CHARACTER* 6 VRN AME ( 6 ) , GENAME ( 6 ) , THI S 
LOGICAL* 1 FINISH ,COMBND,POPULN 

INTEGER OMEGA , VAL , FOOD , TERMOT , TERMIN , AN SR , LAST , DEFALT ,U SER , SAVE 
COMMON/REAL 1/DEFVAL , NUMBER, AMOUNT , GENVAR ,USRVAL ,TR, MCONST , 

*TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA ,MAS ,MFP ,MH ,WWJ ,EBJ ,PJ ,WBJ , 
*ABJ ,REJ ,BOJ ,THJ ,RWJ ,MBJ ,TOTMBJ ,TOTMFP ,MTRJ ,MTP ,MWT ,MPS ,NP ,TBE 
COMMON/CHAR 2/FDNAME , WHICH , LINE ,VRNAME , GENAME , THIS 
COMMON/LOG4/FINISH ,COMBND ,POPULN 

COMMON/INT 3/OMEGA , VAL , FOOD , TERMOT , TERMIN , ANSR , L AST , DEFALT ,U SER 
C 

LAST =6 
SAVE =12 
OMEGA =11 
USER =2 
DEFALT = 1 
TERMOT=6 
TERMIN=5 
FINISH = .FALSE. 

COMBND = .FALSE. 

POPULN = .FALSE. 

C 

WRITE (TERMOT, 5) 

5 FORMAT (' * , 1 THIS PROGRAM IS BASED 

X ON THE PRELIMINARY SCENARIO',/, 

X' ANALYSIS MODEL BY ROBERT M. MASON',//, 

X' TO OBTAIN A DESCRIPTION OF A VARIABLE, TYPE IN— ? VARIABLE ',/ , 
X' NAME— WHEN INPUT IS ASKED FOR') 

OPEN (UNIT=10,NAME=’ MASON .DAT' ,TYPE='OLD' ) 

9 WRITE (TERMOT, 10) 

10 FORMAT (’ ' , ' DO YOU WISH TO RUN A COMBINATION OF FOODS FOR',/, 

* ' BREAKEVEN ANALYSIS , 1=YES , 0=NO ' ) 

C 

READ (TERMIN, 15) LINE 
15 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ. ' ? ' ) THEN 
CALL DESCRT (LINE, 9) 

ELSE 

MODE =1 _ 

CALL MFREIN (LINE , MODE, ANSR ,NVAR, 9) 

END IF “ 

C 

IF (ANSR.EQ. 1) THEN 
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C RUN COMBINATION 

COMBND = .TRUE. 

OPEN (UNIT=12,NAME=' COMB. DAT’ ,TYPE='NEW') 

END IF 
C 

20 CALL PRODUC 

IF A COMBINATION IS RUN, BREAKEVEN CALCULATIONS ARE MANDATORY 
IF (COMBND) THEN 
ANSR = 1 
GO TO 167 
END I F 
C 

159 WRITE (TERMOT ,160) 

160 FORMAT (' 1 , ' DO YOU WISH BREAKEVEN CALCULATIONS , 1=YES , 2=NO ' ) 
READ (TERMIN ,165) LINE 

165 FORMAT (A80) 

C 

IF (LINE(1:1) .EQ. •?• ) THEN 
CALL DESCRT (LINE, 159) 

ELSE 

MODE = 1 

CALL MFREIN (LINE , MODE , ANSR ,NVAR ,159 ) 

END IF 
C 

167 IF (ANSR.EQ. 1) THEN. 

WRITE (TERMOT, 168) 

168 FORMAT (’ BREAKEVEN CALCULATIONS ' ) 

IF ( . NOT . POPULN ) THEN 

169 WRITE (TERMOT, 170) 

170 FORMAT (’ ENTER POPULATION 1 ) 

READ (TERMIN, 172) LINE 

172 FORMAT (A80) 

C 

IF (LINE ( 1: 1) .EQ. ' ? 1 ) THEN 
CALL DESCRT (LINE, 169) 

ELSE 

MODE = 0 

CALL MFREIN (LINE , MODE ,NP ,NVAR ,169 ) 

END IF 
END IF 
C 

CALL BREAK 

WRITE (SAVE, 200) FDNAME (FOOD) ,MTRJ ,MBJ ,MBD , ABJ 
200 FORMAT (All ,F6 . 2 ,F7 . 3 ,F7 . 3 ,F6 . 2) 

END IF 
C 

IF (.NOT .COMBND) THEN 
C 

599 WRITE (TERMOT, 6 00) 

600 FORMAT (' ' , ' DO YOU WISH TO RUN THE PROGRAM AGAIN , 1=YES , 0=NO • ) 

READ (TERMIN, 602) LINE 

602 FORMAT (A80) 

C 

IF (LINE (Is 1) .EQ. ' ?' ) THEN 
CALL DESCRT (LINE, 599) 

ELSE 
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MODE = 1 

CALL MFREIN (LINE , MODE ,ANSR ,NVAR, 599) 

END IF 
C 

IF (ANSR.EQ.l) THEN 
GO TO 20 
END IF 
ELSE 

POPULN = .TRUE. 

605 WRITE (TERMOT , 610 ) 

610 FORMAT (’ ' ANOTHER FOOD? , 1=YES , 0=NO • ) 

READ (TERMOT, 6 15) LINE 
615 FORMAT (A80) 

C 

IF (LINE (1: 1) .EQ. 1 ? 1 ) THEN 
CALL DESCRT (LINE, 605) 

ELSE 

MODE = 1 

CALL MFREIN (LINE , MODE ,ANSR,NVAR, 605) 

END IF “ 

C 

625 IF (ANSR.EQ.0) THEN 

FINISH = .TRUE. 

ENDFILE SAVE 
CLOSE (UNIT=12) 

OPEN (UNIT=12 ,NAME= 1 COMB .DAT' ,TYPE='OLD’) 

ALLMBD = 0 

ALLMTR = 0 

ALLMBJ = 0 

ALLABJ = 0 

690 READ (SAVE ,700 ,END=710 ) MTRJ ,MBJ ,MBD , ABJ 

700 FORMAT ( 11X,F6 . 2 ,F7 . 3 ,F7 . 3 ,F6 . 2) 

ALLMTR = ALLMTR + MTRJ 

ALLMBJ = ALLMBJ + MBJ 

ALLMBD = ALLMBD + MBD 

ALLABJ = ALLABJ + ABJ 

GO TO 690 
C 

710 FDNAME (FOOD) = 'DIET ' 

MBD = ALLMBD 
MTRJ = ALLMTR 
MBJ = ALLMBJ 
ABJ = ALLABJ 
C 

WRITE (TERMOT, 7 20) 

720 FORMAT (' BREAKEVEN CALCULATIONS FOR DIET') 
CALL BREAK 
C 

ELSE 

GO TO 20 
END IF 
END IF 

CLOSE (UNIT=10 ) 

CLOSE (UNIT=12) 

STOP 

END 

C 
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BLOCK DATA 

REAL DEFVAL (6,11) , NUMBER , AMOUNT , GENVAR (6,2) ,USRVAL (6 ,11) ,TR, 

X MCONST ,TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA,MAS ,MFP,MH , 

XWW J , EB J , P J , WB J , AB J , REJ , BO J , TH J , RW J , MB J , TOTMBJ , TOTMFP , 

X MTRJ ,MTP ,MWT ,MPS ,NP,TBE 
CHARACTER*13 FDNAME ( 11) , WHICH 
CHARACTER* 80 LINE 

CHARACTER* 6 VRNAME ( 6 ) , GENAME ( 6 ) , THI S 
LOGICAL* 1 FINISH ,COMBND , POPULN 

INTEGER OMEGA, VAL , FOOD , TERMOT , TERMIN , ANSR, LAST , DEFALT , USER 
COMMON/REAL 1/DEFVAL , NUMBER , AMOUNT , GENVAR ,USRVAL ,TR, MCONST , 

*TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA ,MAS ,MFP ,MH ,WWJ ,EBJ ,PJ ,WBJ r 
*ABJ , REJ ,BOJ f THJ ,RWJ ,MBJ , TOTMBJ f TOTMFP ,MTRJ ,MTP ,MWT ,MPS r NP ,TBE 
COMMON/CHAR2/FDNAME , WHICH , LINE , VRNAME ,GENAME f THIS 
COMMON/LOG 4/FI NISH ,COMBND , POPULN 

COMMON/INT 3/OMEGA , VAL , FOOD , TERMOT , TERMIN , ANSR , L AST , DEFALT , U SER 
DATA FDNAME/’ DRY BEANS ' f ' PEANUT BUTTER' CABBAGE 
X' CARROTS TOMATOES POTATOES ' GREEN BEANS 

X' LETTUCE MELONS ' , ' PEAS WHEAT '/ 

DATA GENVAR/ 1., 0,0,15. 88,83. 2, 122. ,0,0,0, 0,0,0/ 

DATA GENAME/’TR ’,'MFP ’,'MH ','MPSC ', 

X'MA ' , 'MAS '/ 

DATA VRNAME/ ' TH ' , ' TBD ' , ' TBF ' , ' EBD ' , 

X ' EBF ', 'MCONST'/ 

DATA DEFVAL/ 47.0,43.1,204.3,21.1,24.3, .217, 

X110. 0,49. 5, 355. 0,8. 2, 9. 4,. 139, 30. 0,10. 4, 180. 8, 9. 9, 172.1, .059, 
X80. 1,43. 6, 315. 6, 21. 3, 154. 2, .138,215,7.9,113.5,6.2,95.0, .073, 
X120. 0,20. 2, 134. 6, 13. 7, 97. 9,. 150, 60, 108. 9, 920. 6, 26. 3, 305. 2,. 113, 
X28. 0,11. 6, 221. 1,8. 5, 161. 3, .052,107.0,39.9,396.5,19.9,298.9, .083, 
X50. 0,15. 6, 99. 9,. 6, 3. 7, .156, 196. 0,148. 4, 505. 3, 58. 5, 67. 2,. 294/ 

END 


SUBROUTINE PRODUC 

THIS SUBROUTINE PERFORMS CULTIVAR CALCULATIONS 

REAL DEFVAL (6,11) , NUMBER , AMOUNT , GEN VAR ( 6 , 2) ,USRVAL (6 ,11) , TR, 

X MCONST , TBD J ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA,MAS ,MFP ,MH , 

XWWJ , EB J , P J , WB J , AB J , REJ , BO J , THJ , RWJ ,MB J , TOTMBJ , TOTMFP , 

X MTRJ ,MTP ,MWT ,MPS ,NP ,TBE 
CHARACTER* 13 FDNAME ( 11) , WHICH 
CHARACTER* 80 LINE 

CHARACTER* 6 VRNAME ( 6 ) , GENAME ( 6 ) , THI S 
LOGICAL* 1 FINISH, COMBND, POPULN 

INTEGER OMEGA, VAL ,FOOD , TERMOT , TERMIN , ANSR , LAST , DEFALT ,USER 
COMMON/REAL 1/DEFVAL , NUMBER, AMOUNT, GENVAR ,USRVAL ,TR,MCONST , 

*TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA ,MAS ,MFP ,MH ,WWJ ,EBJ ,PJ ,WBJ , 

*AB J , REJ ,BOJ ,THJ ,RWJ ,MBJ , TOTMBJ , TOTMFP ,MTRJ ,MTP ,MWT ,MPS ,NP ,TBE 
COMMON/CHAR 2/FDNAME , WHICH , LINE , VRNAME , GENAME , THI S 
COMMON/LOG 4/FINISH , COMBND , POPULN 

COMMON/INT 3/OMEGA , VAL , FOOD , TERMOT , TERMIN , ANSR , LAST ,DEFALT ,USER 
C 

9 WRITE (TERMOT, 10) 

10 FORMAT (' ' , 'WOULD YOU LIKE A LIST OF FOODS STORED IN THE LIBRARY, 
X1=YES ,0=NO' ) 

READ (TERMIN , 20 ) LINE 
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o o o on 


20 


FORMAT (A80) 

IF (LINE (Is 1) .EQ. 1 ? ' ) THEN 
CALL DESCRT (LINE, 9) 

ELSE 

MODE = 1 
NVAR = 1 

CALL MFREIN (LINE, MODE, ANSR, NVAR,9) 

END IF 
C 

IF (ANSR.EQ. I) THEN 

WRITE (TERMOT ,30) (FDNAME(I) ,1=1, OMEGA) 

30 FORMAT (' ',A13) 

END IF 
C 

49 WRITE (TERMOT, 50) 

50 FORMAT (' ENTER A SINGLE FOOD NAME') 

READ (TERMIN , 55) LINE 

55 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ.'?') THEN 
CALL DESCRT (LINE, 49) 

ELSE 

WHICH = LINE (1:13) 

END IF 
C 

DO 57 1 = 1, OMEGA 

IF (WHICH .EQ .FDNAME (I ) ) THEN 
FOOD=I 
GO TO 35 
ENDIF 

57 CONTINUE 

C 

WRITE (TERMOT, 60) 

60 FORMAT (' ','FOOD NOT IN LIBRARY',/) 

GO TO 9 
C 

35 WRITE (TERMOT, 40) WHICH 

40 FORMAT (' ', 'WOULD YOU LIKE A LIST OF THE VARIABLES AND THEIR 

XDEFAULT VALUES FOR',/,' ' , A13 , ' IN THE LIBRARY , ENTER 1 FOR YES, 
X0 FOR NO' ) 

READ (TERMIN, 45) LINE 
45 FORMAT (A80) 

IF ? THEN GETS DESCRIPTION OF VARIABLE 
IF (LINE (1: 1) .EQ. '? ' ) THEN 
CALL DESCRT (LINE ,35) 

ELSE 

MODE = 1 

NVAR = 1 
RETURNS ANSR 

CALL MFREIN (LINE , MODE , ANSR , NVAR , 35 ) 

END IF 


IF (ANSR.EQ. 1) THEN 
C 

70 DO 74 1=1, LAST 
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WRITE (TERMOT , 71) VRNAME (I ) , DEFVAL ( I , FOOD ) 

71 FORMAT (' ' ,A6,4X,F5.1) 

74 CONTINUE 
C 

END IF 
C 

DO 75 VAL = 1 ,LAST 

USRVAL (VAL, FOOD) = DEFVAL (VAL , FOOD ) 

75 CONTINUE 
C 

77 WRITE (TERMOT ,80) 

80 FORMAT (' 'WOULD YOU LIKE TO CHANGE A DEFAULT VALUE , 1=YES , 0 

X=NO 1 ) 

READ (TERMIN ,85) LINE 
85 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ. ' ? ' ) THEN 
CALL DESCRT (LINE, 77) 

ELSE 

MODE = 1 
NVAR = 1 

CALL MFREIN (LINE , MODE ,ANSR , NVAR, 77 ) 

END IF “ 

C 

C 

IF (ANSR.EQ. 1) THEN 

89 WRITE (TERMOT, 90) 

90 FORMAT (' ', 'WHICH VARIABLE?') 

READ (TERMIN, 95) LINE 

95 FORMAT (A80) 

C 

IF (LINE(1: 1) .EQ. '?' ) THEN 
CALL DESCRT (LINE, 89) 

ELSE 

THIS = LINE (1:6) 

END IF 
C 

98 WRITE (TERMOT, 100) 

100 FORMAT (' ',' ENTER VALUE ' ) 

READ (TERMIN, 101) LINE 

101 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ. ' ? ' ) THEN 
CALL DESCRT (LINE, 98) 

* ELSE 

MODE = 0 

CALL MFREIN (LINE , MODE , NUMBER, NVAR, 98) 

END IF 
C 

DO 110 1=1, OMEGA 

IF (THIS. EQ. VRNAME (I) ) THEN 
USRVAL (I , FOOD ) =NUMBER 
GO TO 77 
END IF 
C 

110 CONTINUE 
C 
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C ERROR ROUTINE 

WRITE (TERMOT,115) 

115 FORMAT (' ' , 'VARIABLE NOT IN FILE') 

GO TO 35 
ENDIF 
C 

C VARIABLES SET TO DEFAULT OR USER DEFINED VALUES 
TH J =U SRVAL ( 1 , FOOD ) 

TBD J =U SRVAL ( 2 , FOOD ) 

TBFJ =U SRVAL ( 3 , FOOD ) 
y EBDJ =U SRVAL (4 ,FOOD) . 

EBFJ =U SRVAL ( 5 , FOOD ) 

MCONST = U SRVAL (6, FOOD) 

C 

* , 139 WRITE (TERMOT, 140) ,FDNAME (FOOD) 

140 FORMAT (' ' , ' TO RUN PROGRAM ENTER DAILY REQUIRMENT OF ',A13,/, 
X' IN GRAMS') 

READ (TERMIN , 141) LINE 

141 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ.'?') THEN 
CALL DESCRT (LINE, 139) 

ELSE 

MODE = 0 

CALL MFREIN (LINE, MODE, AMOUNT, NVAR, 139) 

ENDIF 

C 

C EQUASIONS FOR CALCULATIONS 
REJ = AMOUNT 
WEJ = 1 - EBDJ / EBFJ 

WWJ = 1 - ( (TBDJ - EBDJ) / (TBFJ - EBFJ)) 

EBJ = EBFJ / TBFJ 
PJ = TBFJ 

WBJ = WEJ * EBJ + WWJ * ( 1-EBJ) 

ABJ = REJ / (PJ * EBJ) 

BOJ = REJ / EBJ 

MBJ = (BOJ * THJ / 2) / 1000. 

RWJ = BOJ * (1-EBJ) 

HARWWJ = RWJ * WWJ 
C 

WRITE (TERMOT ,145) 

145 FORMAT (' ' FOOD 10X REJ 5X TH • , 

X 4X , ' TBDJ ' , 3X , ' TBFJ ' , 3X , ' EBDJ 1 , 

1 X3X, 'EBFJ' , 3X , ' EBJ ' ,3X, 'BOJ’ , 5X, 'ABJ' ,3X, 'WEJ' ) 

WRITE (TERMOT, 150) FDNAME (FOOD) , REJ , THJ , TBDJ , 

X TBFJ ,EBDJ ,EBFJ ,EBJ ,BOJ 
X ABJ WEJ 

150 FORMAT (' ' , A13 ,F5 . 1 , 2X ,F5 . 0 , 2X , 4 (F5 . 1 , 2X) , F4 . 3 , IX ,F7 . 1 , 
X1X,F6.2,1X,F4.3) 

C 

RETURN 

END 

C 

C 

C 

SUBROUTINE BREAK 

C THIS SUBROUTINE PERFORMS BREAKEVEN CALCULATIONS 
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c 

REAL DEFVAL (6,11) , NUMBER , AMOUNT ,GEN VAR ( 6 , 2) ,USRVAL ( 6 , 11) ,TR, 

X MCONST ,TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA,MAS ,MFP ,MH , 

XWW J , EB J , P J , WB J , AB J , RE J , BO J ,THJ , RW J ,MBJ , TOTMBJ , TOTMFP , 

X MTRJ ,MTP ,MWT ,MPS ,NP ,TBE 
CHARACTER* 13 FDNAME (11) , WHICH 
CHARACTER* 80 LINE 

CH ARACTER* 6 VRN AME ( 6 ) , GENAME ( 6 ) , THI S 
LOGICAL* 1 FINISH ,COMBND ,POPULN 

INTEGER OMEGA, VAL , FOOD , TERMOT , TERMIN , ANSR,LAST , DEFALT ,USER 
COMMON/REAL 1/DEFVAL , NUMBER, AMOUNT ,GENVAR,USRVAL ,TR , MCONST , 

*TBDJ ,TBFJ ,EBDJ ,EBFJ ,WEJ ,MBD ,MPSC ,MA,MAS ,MFP ,MH ,WWJ ,EBJ ,PJ ,WBJ , 
*ABJ ,REJ ,BOJ ,THJ ,RWJ ,MBJ , TOTMBJ , TOTMFP ,MTRJ ,MTP ,MWT ,MPS ,NP ,TBE 
COMMON/CHAR2/FDNAME , WHICH , LINE , VRNAME , GENAME , THIS 
COMMON/LOG4/FINISH ,COMBND ,POPULN 

COMMON/INT 3/OMEGA ,VAL ,FOOD , TERMOT , TERMIN , AN SR, LAST , DEFALT ,USER 
C 

179 WRITE (TERMOT, 180) 

180 FORMAT (' • , ' DO YOU WISH A LIST OF VARIABLES AND DEFAULT VALUES, 
XO=NO , 1=YES 1 ) 

READ (TERMIN, 18 5) LINE 
185 FORMAT (A80) 

C 

IF (LINE (Is 1) .EQ. •? • ) THEN 
CALL DESCRT (LINE, 179) 

ELSE 

MODE =1 

CALL MFREIN (LINE , MODE ,ANSR ,NVAR, 179 ) 

END IF 
C 

IF (ANSR.EQ. 1) THEN 
DO 210 1=1, LAST 

WRITE (TERMOT, 195) GENAME (I ) ,GENVAR (I , DEFALT) 

195 FORMAT (' ' , A6 , 3X ,F6 . 2) 

210 CONTINUE 

ENDIF 
C 

C SETS DEFAULT VALUES TO USER VALUES 
DO 220 1=1, LAST 

GENVAR (I ,USER) =GENVAR (I , DEFALT) 

220 CONTINUE 
C 

225 WRITE (TERMOT, 2 30) 

230 FORMAT (' ' ,'DO YOU WISH TO CHANGE THE VALUE OF A VARIABLE, 
X1=YES , 0=NO' ) 

READ (TERMOT, 240) LINE 
240 FORMAT (A80) 

C 

IF (LINE (1:1) .EQ. ’?* ) THEN 
CALL DESCRT (LINE, 225) 

ELSE 

MODE = 1 

CALL MFREIN (LINE , MODE ,ANSR ,NVAR , 225 ) 

END IF 
C 

IF (ANSR.EQ. 1) THEN 
249 WRITE (TERMOT, 250) 
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o o on on 


250 

FORMAT (' ' WHICH 

VARIABLE ' ) 


READ (TERMIN ,260) 

LINE 

260 

r* 

FORMAT (A80) ■ 



IF (LINE (1:1) .EQ. '?' 

) THEN 


CALL DESCRT (LINE, 24 9) 
ELSE 


THIS = LINE (1:6) 
END IF 



c 


•0 

269 

WRITE (TERMOT, 270) 


270 

FORMAT ( ' ' , ’ENTER VALUE 



READ (TERMIN, 27 2) LINE 


272 

FORMAT (A80) 


C 



IF (LINE (1:1) .EQ.'?') THEN 
CALL DESCRT (LINE, 269) 

ELSE 

MODE = 0 

CALL MFREIN (LINE , MODE , NUMBER ,NVAR ,269 ) 

END IF 
C 

DO 280 1=1, LAST 

IF (THIS .EQ.GENAME (I) ) THEN 
GENVAR (I ,USER) =NUMBER 
GO TO 225 
ENDIF 

280 CONTINUE 

C 

WRITE (TERMOT ,285) 

285 FORMAT (' VARIABLE NOT IN FILE’) 

GO TO 179 
ENDIF 

SETS VARIABLES EQUAL TO USER VALUES 
TR = GENVAR (1, USER) 

MFP = GENVAR (2, USER) 

MH = GENVAR (3, USER) 

MPSC = GENVAR (4, USER) 

MA = GENVAR (5, USER) 

MAS = GENVAR (6, USER) 

IF THE PLANT COMBINATION IS NOT COMPLETED 
IF (.NOT. FINISH) THEN 
MBD = MBJ * MCONST 

MTRJ = (EBDJ * 365.25 * 1.53 * ABJ) /1000. 

END IF 

EQUASIONS FOR BREAKEVEN CALCULATIONS 

MPS = (MPSC) + (MAS * (ABJ * NP) ** .6) + (MA * ABJ * NP) 
MWT = 10. * MBD * NP * TR 
TOTMFP = MFP * (NP /10.) ** .415 
TOTMB = NP * MBJ 

MTP = TOTMB + MH + TOTMFP + MPS + MWT 
TBE = MTP / (MTRJ * NP) 

C 

WRITE (TERMOT, 290) 
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290 FORMAT {' • , 3X , ' FOOD* , 28X, 'MTRJ' , 5X, 'MB' , 8X, 'MPS ', 8X , 'MWT* , 6X, 
X'TBE',/) 

WRITE (TERMOT, 200) FDNAME (FOOD) , MTRJ ,MBJ ,MPS , MWT ,TBE 
200 FORMAT (' ' , 3X , A13, 18X , F6. 2 , 2X , F7 . 3 , 2X,F8 . 1 ,2X ,F8 . 1 f 2X ,F7 . 1 ,/) 
C 

RETURN 

END 

C 

SUBROUTINE DESCRT (LINE , * ) 

C THIS SUBROUTINE RETURNS THE DESCRIPTION OF ANY VARIABLE WHEN A '?' 
C VARIABLE NAME IS ENTERED 
C 

CHARACTER* 80 LINE 
CHARACTER* 90 RECORD 
CHARACTER* 6 VARBLE 
CHARACTER* 83 SCRIPT 
INTEGER FILE, TERMOT 
C 

FILE = 10 
TERMOT =6 
C 

C READS THE FILE THAT CONTAINS DESCRIPTIONS 
5 READ (FILE, 10 ,END=25) VARBLE , SCRIPT 

10 FORMAT (A6 , IX , A83) 

C 

C COMPARES WITH FILE TO FIND DESCRIPTION 
IF (LINE (2: 7) .EQ. VARBLE) THEN 
WRITE (TERMOT, 15) VARBLE , SCRIPT 
15 FORMAT (' ' ,A6 , 3X ,A83) 

GO TO 40 
C 

C READS ANOTHER 
ELSE 

GO TO 5 
C 

END IF 
C 

C ERROR ROUTINE 
25 WRITE (TERMOT, 30) 

30 FORMAT (’ DESCRIPTION OF VARIABLE NOT FOUND') 

C 

C SETS POINTER AT BEGINNING OF FILE 
40 REWIND FILE 

RETURN 1 
END 
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